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Abstract: We give explicit formulae for the 0(e) and 0(e 2 ) contributions to the unrenor- 
malised three loop QCD corrections to quark and gluon form factors. These contributions 
have at most transcendentality weight eight. The 0(e) terms of the three-loop form fac- 
tors are required for the extraction of the four-loop quark and gluon collinear anomalous 
dimensions. The 0(e 2 ) terms represent an irreducible contribution to the finite part of the 
form factors at four- loops. For the sake of completeness, we also give the contributions to 
the one and two loop form factors to the same transcendentality weight eight. 
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The form factors are fundamental ingredients for many precision calculations in QCD. 
These basic building blocks describe the coupling of an external, colour-neutral off-shell 
particle to a pair of partons: the quark form factor is the coupling of a virtual photon to 
a quark-antiquark pair, while the gluon form factor is the coupling of a Higgs boson to a 
pair of gluons through an effective Lagrangian. 

The form factors are phenomenologically important and appear directly as virtual 
higher-order corrections in coefficient functions for the inclusive Drell-Yan process [1-3] 
and the inclusive Higgs production cross section [3-6]. The form factors also display a 
non-trivial infrared pole structure which is determined by the infrared factorisation for- 
mula. This implies that their infrared pole coefficients can be used to extract fundamental 
constants such as the cusp anomalous dimensions which control the structure of soft di- 
vergences and the collinear quark and gluon anomalous dimensions. In fact, the cusp 
anomalous dimensions were first obtained to three loops from the asymptotic behaviour 
of splitting functions [7, 8]. However, it was the calculation [9, 10] of the pole terms of 
the three-loop form factors (and finite plus subleading terms in the two-loop and one-loop 
form factors [11-13]), which led to the derivation of the three-loop collinear anomalous 
dimensions [9, 14, 15]. 

The infrared factorisation formula for a given form factor (or more generally for a given 
multi-leg amplitude) at a certain number of loops involves infrared singularity operators 
acting on the form factor evaluated with a lower number of loops. These infrared singularity 
operators contain explicit infrared poles 1/e 2 and 1/e. Therefore, the computation of the 
finite contribution to any re-loop form factor relies on contributions from (n — m)-loops 
evaluated to C(e 2m ). 

At present, the state of the art is at the three- loop level for the massless quark and 
gluon form factors. There are 22 master integrals shown in Fig. [I], of which 14 are genuine 
three-loop vertex functions (A^-type) , 4 are three- loop propagator integrals (B^-type) and 
4 are products of one-loop and two-loop integrals (C^-type). In this notation, the index t 
denotes the number of propagators, and i is simply enumerating the topologically different 
integrals with the same number of propagators. Expressions for the form factors in terms of 
the 22 independent master integrals, and valid for any value of the dimension D, are given in 
Ref. [16]. The -B^-type integrals were computed to finite order in [17, 18] and supplemented 
by the higher order terms in [19]. Explicit expansions of the j4^j-type integrals were 
obtained in Refs. [20-23] using Mellin-Barnes techniques. They enabled the evaluation of 
the three-loop form factors up to and including the finite contributions [16, 23, 24]. The 
deepest pole contribution is of 0(l/e 6 ). Correspondingly, the finite terms are of at most 
transcendent ality weight six, that is terms such as 7r 6 (£f) or £|. 

More recently [25] , 20 of the three- loop master integrals have been re-evaluated up to 
transcendent ality weight eight using dimensional recurrence relations [26, 27] and analytic 
properties of Feynman integrals (the DRA method [28]). Expressions for the two remaining 
integrals, i?s,i and Cs,i, can be obtained from Refs. [28] and [13] respectively. Once the 
same normalisation and basis set of multiple zeta values is used, Ref. [25] confirms the 
earlier result of Ref. [29] for Aq^- On the other hand, we confirm a certain subset of mas- 
ter integrals (-66,2, B$ t i, Aj^, A^^, A%j_, -Ag.i, ^9,2, ^9,4) from [25, 28] up to coefficients 
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Figure 1: Master integrals for the three-loop form factors. Labels in brackets indicate the naming 
convention of Ref. [25] . 

corresponding to weight eight numerically to a precision of one per-mille or better using 
MB.m [30] and FIESTA [31, 32]. All other of the 22 master integrals we even confirm ana- 
lytically through to weight eight by expanding the closed form in terms of hypergeometric 
functions given in [20, 21] using the HypExp package [33]. 
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All master integrals are therefore known up to transcendentality weight eight i.e. terms 
including 7r 8 (Q), C2C3 > C3C5 as wen as t ne multiple zeta value Cs,3 (or equivalently C-6,-2)- 
The multiple zeta values are defined by (see e.g. [34] and references therein) 

00 i\— 1 ik-i~ 1 k / \i ■ 

«i = EE- £ <u 

ii=li 2 =l «fe=l j=l z j 

Specifically, C-6,-2 is related to Cs,3 by [34, 35] 

9 3 781 8 

C-6,-2 - ^C5,3 - ^CsCs + • (2) 

The numerical values of the transcendental constants up to weight eight are: 

C 3 = 1.2020569031595942854 . . . , Cs = 1.0369277551433699263 . . . , 
C 7 = 1.0083492773819228268 . . . , Cs,3 = 0.037707672984847544011 .... 

The new results on the higher order terms in the master integrals enable the computation 
of the three-loop form factors through to 0(e 2 ) which is an intrinsic component for the 
four-loop evaluation of the form factors. This is the topic of this Letter and we give 
explicit formulae for the 0(e) and C(e 2 ) contributions to the unrenormalised three loop 
form factors. 

The form factors are the basic vertex functions of an external off-shell current (with 
virtuality q 2 = S12) coupling to a pair of partons with on-shell momenta p\ and p2- One 
distinguishes time-like (s±2 > 0, i.e. with partons both either in the initial or in the final 
state) and space-like (si2 < 0, i.e. with one parton in the initial and one in the final state) 
configurations. The form factors are described in terms of scalar functions by contracting 
the respective vertex functions (evaluated in dimensional regularization with D = 4 — 2e 
dimensions) with projectors. For massless partons, the full vertex function is described 
with only a single form factor. 

The quark form factor is obtained from the photon-quark-antiquark vertex by 

f q = -^rhw Tr & r «M, (3) 

while the gluon form factor relates to the effective Higgs-gluon-gluon vertex Fgg as 

T g _ Pi • P2 9^ ~ Pl,yP2,v ~ Pl,»P2,y ^ , , 

2(1 -e) 99 • { > 

The form factors are expanded in perturbative QCD in powers of the coupling constant, 
with each power corresponding to a virtual loop. We denote the unrenormalized form 
factors by T a and the renormalized form factors by F a with a = q,g. 

At tree level, the Higgs boson does not couple either to the gluon or to massless quarks. 
In higher orders in perturbation theory, heavy quark loops introduce a coupling between 
the Higgs boson and gluons. In the limit of infinitely massive quarks, these loops give rise 
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to an effective Lagrangian [36] mediating the coupling between the scalar Higgs field and 
the gluon field strength tensor: 



Ant - --HF£ u F attlu 



(5) 



The coupling A has inverse mass dimension. It can be computed by matching [37, 38] the 
effective theory to the full standard model cross sections [5]. 

Direct evaluation of the Feynman diagrams at the appropriate loop order yields the 
bare (unrenormalised) form factors, 



J^(a b s ,s 12 ) 




on -pq 
D e J- n , 



(6) 
(7) 



where //q is the mass parameter introduced in dimensional regularisation to maintain a 
dimensionless coupling in the bare Lagrangian density and where 



S e 
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(47T) 



with the Euler constant 7 = 0.5772 . 



(8) 



The one-loop and two-loop form factors were computed in many places in the litera- 
ture [9-13]. All-order expressions in terms of one-loop and two- loop master integrals are 
given in [13]. Explicit expressions for the one- and two- loop form factors through to 0(e 5 ) 
and 0(e 3 ) respectively are given already in [16]. To determine the finite piece at the 
four-loop level, these form factors are needed to one higher power in e, and for the sake of 
completeness, we quote them here. At one-loop, 
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and at two-loops 
J-S = I 1 + . . . + J-% 



2 \e J 



+ C 2 F 



+ - 



637631 
128 



528C 5 ,3 + 



27204C 7 34001(5 
7 10 



481913C 3 33248CsC5 



24 



15 



- 4 - 



36359C; 2 



3 , 95559C2 , , 2257C 2 C.3 4576C 2 C| 248023C 2 2 
H — iyoC2<,5 H ^ p; 



32 



80 



| 5109C 2 2 C 3 55623CI 653901(1 



5 ' 140 

11630115085 
839808 



700 



+ 264C 5 3 - 1198 ° Cr + 1214029 ^> + 84520897(3 



21 



270 



8266C 3 C 5 229042CI 58499773(2 829C 2 C 5 94931C 2 Cs , 3029C 2 C 3 2 



81 



23328 



15 



162 



5832 
+ 



14915741(1 66379C|Cs 4843C| 75242C| 



C F N F 



6480 90 30 

996726245 2186C 7 42713C 5 



+ 



419904 



21 



135 



175 

1951625C3 
2916 



+ 



4732C| , 2877653C 2 



81 



+ 



11664 



242C2C5 4589C 2 C3 309181C 2 2 533CIC3 127C| 



15 



81 



3240 



45 



T 9 - 

x 2 — 



+ cl 



1371828689 56155C 7 161266C 5 5108944C 3 I69OC3C5 
^ ; 7. 264C5,3 H — t^z 1- 



209952 



42 



135 



729 



85559Cf 219275C2 IOOIC2C5 H858C2C3 1547C 2 Cl 187733(1 



11 



1944 



+ 



27 



180 



+ 22781C|C3 123079CI 50419(1 



90 



1260 



100 



+ C A N F 



232282297 229( 7 _ 24518C 5 _ 301886C 3 22060C 3 2 



104976 



21 



135 



729 



81 



+ 98791C2 + 342C2C5 + 2978C2C3 _ 40148C 2 2 + 517C 2 2 C 3 + 2167C 2 3 



+ C F N F 



972 5 
19296691 



27 



405 



630 



7776 



45 



81 



648 



-32C2C5 - 



5716C2C3 , 585929C 2 2 6724C 2 2 Cs 2024C 2 3 



27 



+ 



1620 



45 



105 



(10) 



The unrenormalised three- loop quark form factor J| through to (and including) O(e ) 
is given in eq. (5.4) of Ref. [16]. The pole contributions of are also given in eq. (3.7) 
of ref. [9] while the finite parts of the Np, CaN f and C F N F contributions are given in 
eq. (6) of ref. [10]. The finite Npy contribution could already be inferred from [39]. The 
remaining finite contributions are also given in eqs. (8) and (9) of ref. [24]. The 0(e 1 ) and 
C(e 2 ) contributions are given by, 
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Note that last colour factor is generated by graphs where the virtual gauge boson does 
not couple directly to the final-state quarks. This contribution is denoted by Npy and 
is proportional to the charge weighted sum of the quark flavours. In the case of purely 
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electromagnetic interactions, we find, 



(12) 



The unrenormalised three-loop gluon form factor through to (and including) O(e ) is 
given in eq. (5.5) of Ref. [16]. The divergent parts are also given in eq. (8) of ref. [10] while 
the finite contributions are given in eq. (10) of ref. [24]. The 0(e l ) and 0(e 2 ) contributions 
for T?, are given by, 
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The renormalised form factors are directly related to the unrenormalised form factors 
and details on how to extract the renormalised form factors to this order are given in 
section 2 of Ref. [16]. 

In this letter, we computed the three-loop quark and gluon form factors through to 
0(e 2 ) in the dimensional regularisation parameter. These contributions are relevant in the 
study of the infrared singularity structure at four loops. In particular, the 0(e) terms 
of the three-loop form factors are required for the extraction of the four-loop quark and 
gluon collinear anomalous dimensions. The 0(e 2 ) terms contribute to the finite part of 
the infrared-subtraction of the form factors at four loops. It is this infrared-subtracted 
finite part which is relevant for the study of the next-to-next-to-next-to-next-to-leading 
(N 4 LO) Drell-Yan and Higgs production processes. In particular, the 0(e 2 ) three-loop 
contributions represent a finite ingredient to these processes at four-loops. 
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